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ABSTRACT 

We examine dust-to-gas ratio as a function of metallicity for dwarf galaxies 
[dwarf irregular galaxies (dlrrs) and blue compact dwarf galaxies (BCDGs)]. 
Using a one-zone model and adopting the instantaneous recycling approximation, 
we prepare a set of basic equations which describes processes of dust formation 
and destruction in a galaxy. Four terms are included for the processes: dust 
formation from heavy elements ejected by stellar mass loss, dust destruction in 
supernova remnants, dust destruction in star-forming regions, and accretion of 
heavy elements onto preexisting dust grains. Solving the equations, we compare 
the result with observational data of nearby dlrrs and BCDGs. The solution 
is consistent with the data within the reasonable ranges of model parameters 
constrained by the previous examinations. This means that the model is 
successful in understanding the dust amount of nearby galaxies. We also show 
that the accretion rate of heavy element onto preexisting dust grains is less 
effective than the condensation of heavy elements. 



Subject headings: dust, extinction — galaxies: evolution — galaxies: ISM - 
galaxies: irregular 



1. INTRODUCTION 

Interstellar dust is composed of heavy elements made and ejected by stars. Dwek & 
Scalo (1980) demonstrated that supernovae are the dominant source for the formation of 
dust grains. They also showed that the dust is destroyed in supernova shocks (see also 
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McKee 1989, Jones et al. 1994, and Jones, Tielens, & Hollenbach 1996). Thus, the dust is 
formed and destroyed in star-forming galaxies. 

Some recent galaxy-evolution models treat the evolution of total dust mass as well 
as that of metal abundance (Wang 1991; Lisenfeld & Ferrara 1998, hereafter LF98; Dwek 
1998, hereafter D98; Hirashita 1999, hereafter H99). The processes of dust formation and 
destruction by supernovae are taken into account in LF98, in order to explain the relation 
between dust-to-gas mass ratio and metallicity of dwarf irregular galaxies (dlrrs) and 
blue compact dwarf galaxies (BCDGs). We can find detailed mathematical formulations 
to calculate the dust mass in any galactic system in D98, which treats the accretion of 
heavy elements onto preexisting dust grains in molecular clouds in addition to the processes 
considered in LF98. The model by D98 is successfully applied to nearby spiral galaxies in 
H99. 

It is suggested in D98 that the accretion process onto preexisting dust grains is not 
effective in dwarf galaxies because of the absence of dense molecular clouds there. If this 
is true, it is worth estimating the ineffectiveness quantitatively, by which we can obtain 
an information about molecular clouds in dwarf galaxies. Direct observations of molecular 
clouds in some dwarf galaxies have been extensively carried out (e.g., Ohta, Sasaki, & Saito 
1988). However, it is generally difficult to observe extragalactic molecular clouds. Thus, 
theoretical constraints on the nature of molecular clouds are indispensable for investigations 
on the star formation processes in extragalactic objects. 

In this paper, we examine the relation between dust-to-gas ratio as a function of 
metallicity by using a set of equations in H99. The model is one-zone (i.e., the spatial 
distribution of physical quantities within a galaxy is not taken into account) and the 
instantaneous recycling approximation is applied. The main purpose of this paper is to 
apply the model to dwarf galaxies (dlrrs and BCDGs). First of all, in the next section, we 
explain the model equation, through which the dust-to-gas ratio as a function of metallicity 
is calculated. The result is compared with observational data of dlrrs and BCDGs in §3. 
We present a summary in the final section. 

2. REVIEW OF MODEL EQUATIONS 

In order to investigate the dust content in galaxies, H99 derived a set of equations 
describing the dust formation and destruction processes, based on the models in LF98 and 
D98. In H99, a simple one-zone model for a galaxy is adopted to extract global properties 
of galaxies. For the model treating the radial distribution of gas, element, and dust in a 



- 3- 



galaxy, see D98. 



The equation set is written as 



(1) 



dt 
dMi 



-Xnl> + Et- XiW, 



(2) 
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(See eqs. [l]-[3] in McKee 1989, eqs. [6]-[8] in LF98 and eqs. [l]-[3] in H99.) Here, M g is 
the mass of gas. The metal is labeled by % {% = 0, C, Si, Mg, Fe, ■ • •), and and M dji 
denote the total mass of the metal i (in gas and dust phases) and the mass of the metal 
i in the dust phase, respectively. The star formation rate is denoted by ip; E is the total 
injection rate of mass from stars; W is the net outflow rate from the galaxy; Xi is the mass 
fraction of the element i (i.e., Xi = Mj/M g ); Ei is the total injection rate of element i from 
stars; /, is the mass fraction of the element i locked up in dust (i.e., fi = M^/M,). The 
meanings of the other parameters in the above equations are as follows: f m> i is the value 
of the dust mass fraction in the injected material, in other words, the dust condensation 
efficiency in the eject a|]; a refers to the efficiency of the dust destruction during star 
formation [a = 1 corresponds to destruction of only the dust incorporated into the star, 
and a > 1 (a < 1) corresponds to a net destruction (formation) in the star formation]; 
r acc is the accretion timescale of the element i onto preexisting dust particles in molecular 
clouds; tsn is the timescale of dust destruction by supernova shocks; 6 accounts for the dust 
content in the outflow (5 = means no dust in the outflow, while 5 = 1 indicates that the 
outflow is as dusty as the interstellar medium). 

We should comment on the parameter a here. Since the protostellar disk forms dust, 
a < 1 is expected. However, as to the circumstellar dust, the timescale of loss of angular 
momentum through the Poynting-Robertson effect is much shorter than the lifetime of stars 
(e.g., Rybicki & Lightman 1979). This means that the formed dust is lost effectively. Thus, 
we reasonably assume that a = 1 hereafter. The formation of planets also contributes to 
the loss of the dust. 

Here, we adopt the same assumption as LF98 and H99; the instantaneous recycling 
approximation (Tinsley 1980): Stars less massive than mi live forever and the others die 



this formalism, we assume that the condensation efficiency in stellar winds is the same as that in 
supernova ejecta. 
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instantaneously. This approximation allows us to write E and E i} respectively, as 

E = Kil>, (4) 
Ei = (KXt + yM, (5) 

where TZ is the returned fraction of the mass that has formed stars which is subsequently 
ejected into the interstellar space, and 3^ is the mass fraction of the element i newly 
produced and ejected by stars.Q TZ and 3^ can be obtained using the following formulae 
(Maeder 1992): 

TZ — / (m — w m )<j)(m)dm, (6) 

J mi 

y% — mpi(m)(j)(m)dm, (7) 

In equation (|^), (f){m) is the initial mass function (IMF), and m u is the upper mass cutoff of 
stellar mass. The IMF is normalized so that the integral of mcf)(m) in the full range of the 
stellar mass becomes 1. Therefore, <p(m) has a dimension of the inverse square of the mass. 
In equation (|7|), w m is the remnant mass (w m = 0.7M Q for m < 4M Q and w m = 1.4M for 
m > 4M ) and Pi(m) is the fraction of mass converted into the element i in a star of mass 
m. 

Using the above parameters TZ and 3^, and assuming that W is proportional to the 
star formation rate (W = wif)), equations (ffl) - ® become 



1 dM g 
ip dt 

M g dXi 
ifj dt 

M r dV; 



-l + n-w, (8) 

y» (9) 

f^TlXi + y i )-[a-l + K- /3 acc (l - fc) + f3 SN - w(l - 5)}V h (10) 



ip dt 

where T>i = M&^jM g = fcXi, and /3 acc and @sn are, respectively, defined by 

/3 acc = ^- and /3 SN = ^%. (11) 

We can regard /3 acc and /3sn as constant in time (§6.2 and §8.4 of D98). We note that the 
Galactic value shows /3sn ~ 5 (LF98). This value corresponds to r SN ~ 10 8 yr, which is 
consistent with Jones et al. (1994) and Jones, Tielens, & Hollenbach (1996). The relation 
r acc ^ W2 (D98) leads to /3 acc ~ 2/3 SN ~ 10. 



2 K = R and ^ = y(l - R) for the notation in LF98. 
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Combining equations (J|) and fllPp, we obtain the following differential equation of T>i 
function of X,: 

y*j^ = fin,i{nXi + y i )-[a-i + n- p m + /3 SN - w(i - 5)] v, - t^-L, (12) 

where we used the relation, ft = T>i/Xi. 

Here, we take % = O to compare the result with the data in LF98. For further 
quantification, we need to fix the values of IZ and 3^ for the traced element (i = O). 
The values are given in LF98. The reason why LF98 chose oxygen as the tracer is as 
follows: (i) Most of oxygen is produced in Type II supernovae which are also responsible 
for the shock destruction of dust grains; (ii) oxygen is the main constituent of dust grains. 
The point (i) means that the instantaneous recycling approximation may be a reasonable 
approximation for the investigation of oxygen abundances, since the generation of oxygen 
is a massive-star-weighted phenomenon. In other words, results are insensitive to the value 
of m\. According to LF98, we put mi = 1M and m u = 12OM . We use a power-law form 
of the IMF: (j)(m) oc m~ x . The Salpeter IMF is assumed; i.e., x = 2.35 (Salpeter 1955). 
According to LF98, (11, y ) = (0.79, 1.8 x 1(T 2 ) for the Salpeter IMF. After the numerical 
integration of equation (0) by the Runge-Kutta method, we compare the result with 
observational data of dwarf galaxies in the next section. 



3. APPLICATION TO DWARF GALAXIES 

To compare the solution of equation (|T2|) with observational data of dust-to-gas ratio, 
we make an assumption that total mass of dust is proportional to that of oxygen in the 
dust phase. In other words, 

V = J2^ = CV , (13) 

i 

where C is assumed to be constant for all galaxies. According to Table 2.2 of Whittet 
(1992), C ~ 2.2 (the Galactic value). 

We compare the solution of equation (O) with the data in LF98 (see also the references 
therein). The data sets of nearby dlrrs and BCDGs are presented in Tables 1 and 2 of 
LF98, respectively. The observed dust-to-gas ratio is defined as 

V ohs = M d /M m , (14) 

where M d and M H i are the total masses of dust and H I gas, respectively. The dust mass is 
derived from the luminosity densities at the wavelengths of 60 /zm and 100 /im observed by 



- 6- 



IRAS. The dust mass derived from the far-infrared emission is about an order of magnitude 
smaller than the value found from the analysis of the interstellar extinction (Fig. 2 of 
LF98). The presence of cold or hot dust, emitting beyond 100 /xm and below 60 /im may be 
responsible for the discrepancy (LF98). Thus, we should keep in mind that the dust-to-gas 
ratio adopted here is underestimated in this context. However, since we only take into 
account the H I gas as the gas content and do not consider H 2 gas, the dust-to-gas ratio is 
overestimated (by a factor of ~ 2). To sum up, we should keep in mind the uncertainty of 
an order of magnitude for the dust-to-gas ratio derived from the observation (T> ohs ). 

In the following subsections, we compare V calculated by using equations ( |T2"D and ( pT3|) 
with V obs . We focus on the two processes of dust formation: One is the condensation of 
dust from heavy elements ejected by stars, and the other is the accretion onto preexisting 
dust grains. The efficiency of the former process is denoted by /i nj o and that of the latter 
by /9 a cc- The latter process is not taken into account in LF98. For the dependences of the 
relation between dust-to-gas ratio and metallicity on IMF, see LF98 and H99. 



3.1. Dependence on /; n o 

We present the dependence of the result on the value of / inj o in Figure la, in which 
the solid, dotted, and dashed lines show the T>-X relation for / in = 0.1, 0.05 and 0.01, 
respectively. The other parameters are fixed to a = 1, j3 acc = 2/?sn = 10, 5 = 1 [i.e., 
w(l — 8) = 0]. Figure la shows that the larger the efficiency of production of dust from 
heavy elements is, the larger the dust-to-gas ratio becomes. The data points represent the 
relations between T> ohs and Xq of dlrrs and BCDGs in LF98. The filled and open squares 
show the data points of the dlrrs and the BCDGs, respectively. 

In the limit of X Q — > 0, the solution reduces to 

^o-/in,o^o or V~Cf^ X . (15) 

(See also LF98.) This means that T> scales linearly with /i n .o for the extremely low 
metallicity. Thus, /i n) o can be constrained by low-metal galaxies (see also H99). Equation 
(|15|) means that we can constrain the parameter /i n) o by examining the relation between 
Uo/Xq (the fraction of oxygen in the dust phase) and X G . We show this relation in Figure 
lb. The parameter adopted for each line is the same as Figure la. We also show the data 
points of the same sample as Figure la, assuming C = 2.2 for all the galaxies to convert 
pobs j n ^- Roughly speaking, our model predicts that Vo/Xq is constant in the range 
of the dwarf galaxies. This indicates that the low-metal approximation used to derive 
equation ([15]) is applicable to dwarf galaxies. Thus, we can directly constrain the parameter 
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/in,o by dwarf galaxies. From the data points in Figure lb, we see 0.01 <^ V /X ^0.1 or 
0.01 < /i ni o ^ 0.1. We note that this range is consistent with the analyses in H99. However, 
we should keep in mind the uncertainty of the data as described above in this section. 

We also see from Figure lb that the relation between T>o and Xo becomes nonlinear 
in the relatively high- metal region (logXo > —3). The behavior of this nonlinear region 
depends on /3 acc or /3 SN (§3.2). 

3.2. Dependence on f3 acc 

We here investigate the dependence of the solution on (3 acc , which is proportional 
to the accretion efficiency of heavy elements onto the preexisting dust grains (§2). The 
resulting V-Xq relations for various j3 acc are shown in Figure 2a. We show the cases of 
Pace = 5, 10, 20 (the solid, dotted, and dashed lines, respectively), where the relation 
Pace = 2/?sn is fixed. The other parameters are set to /; n> o = 0.05 (~ the center of the range 
constrained in §3.1), a — 1, and 5 = 1. The value of /3 acc is determined by the lifetime 
of molecular clouds (D98). The value /3 acc ~ 10 corresponds to the accretion timescale of 
~ 10 8 yr (§2). The increase of /3 acc means that the accretion of heavy elements onto dust 
becomes efficient. Thus, for a fixed value for the metallicity, dust-to-gas ratio increases as 
/5 acc increases. 

We also present an extreme case of /5 acc = and Psn = 5 (long-dashed line). In this 
case, the accretion onto preexisting dust grains is neglected. We were not able to reproduce 
the observational data of nearby spiral galaxies without taking into account the accretion 
process (D98; H99). However, the solid and long-dashed lines in Figure 2 show that we 
cannot judge whether the accretion process is efficient or not because of the little difference 
between the results with and without the accretion process. Thus, the basic equations of 
LF98, which do not include the term of the accretion were able to explain the observed 
relation between dust-to-gas ratio and metallicity. We note that the accretion process is 
properly considered in D98. 

The ineffectiveness of the accretion process is understood as follows. Two processes 
are responsible for the formation of dust in equation (|3]); the dust condensation from the 
heavy elements ejected by stars and the accretion of the heavy elements onto preexisting 
dust grains. The former is expressed as f\ n ,iEi and the latter as Md,i(l — fi)/r acc (see eq. 
@). We define the following ratio, Af. 

M d ,(l - ft/Tt 
% f -E- 

J m,i- LJ i 



(16) 
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If Ai < 1, the accretion process is ineffective compared with the dust condensation process. 
We will show that Ai < 1 for the sample dwarf galaxies. 

Using the instantaneous recycling approximation, Ai is estimated as 

Now we put i — O. In Figure 2b, we present the relation between fo = Vo/Xo and Xo- 
The values of parameters for each line is the same as Figure 2a. This figure shows that 
we can consider that 1 — fo ~ 1. Moreover, in the range in which we are interested here, 
Vq ~ fin,oXo (eq. fl5fl), and TZX <C y Q (as long as \ogX Q ^ —2 is satisfied). Thus, A Q 



can be approximated by 



.4 



o 



ftaccXt 



O 



If we put /3 acc = 10, and y = 10~ 2 , we obtain logX G ^ —3 for the condition A Q < 1. This 
is consistent with Figure 2b, since the difference between the solid and the long-dashed 
lines is clear for logXo > —3. Thus, if logXo > —3 is satisfied, the dust accretion process 
is more ineffective than the dust condensation process. Actually, even for logXo ~ —2.5, 
the difference is within the typical error of the observed values. 

The ineffectiveness of the accretion process results from the low metallicity. Thus, in 
galaxies with high metallicity, the accretion process becomes important. Indeed, D98 and 
H99 showed that the process is effective in spiral galaxies, whose metallicity is much larger 
than the dwarf galaxies (H99). 

Since we can reproduce the relation between dust-to-gas ratio and metallicity of dwarf 
galaxies without considering the accretion process, D98 suggested that the accretion process 
is not efficient in dwarf galaxies because of the lack of dense molecular clouds. This may 
be true, but seeing that dlrrs and BCDGs show high star formation efficiency (Sage et al. 
1992; Israel, Bontekoe, & Kester 1996), there may be a large amount of dense molecular gas 
in dwarf irregular galaxies. Indeed, we have shown in Figure 2 that the observed relation 
can be explained even if the efficiency of the accretion process /3 acc is as high as that in the 
spiral galaxies considered in H99. 



4. SUMMARY AND DISCUSSIONS 



Based on the models proposed by LF98 and D98, we have examined the dust content 
in dlrrs and BCDGs. The basic equations which describes the changing rate of dust-to-gas 
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ratio include the terms of dust formation from heavy elements ejected by stars, destruction 
by supernova shocks, destruction in star-forming regions and accretion of elements onto 
preexisting dust grains (§2). This accretion process is important in molecular clouds, where 
gas densities are generally high. The results are compared with the observed values of dlrrs 
and BCDGs. Though the degeneration of the parameter and observational error makes 
it impossible to determine each of the parameter precisely, we were able to constrain the 
parameters to some extent. 

The efficiency of dust production from heavy element (denoted by /i n ,i) can 
be constrained by the galaxies with low metallicity (§3.1). The reasonable range is 
0.01 <^ f ilh i <^ 0.1, which is consistent with H99. Thus, it is possible to understand the dust 
amount in dwarf systems as well as that in spiral systems through the model in this paper. 

As for the nearby spiral galaxies, unless we take into account the accretion process of 
heavy element onto the preexisting dust particles, we cannot explain the observed relation 
between dust-to-gas ratio and metallicity (D98; H99). For the dwarf galaxies, however, we 
can explain the data without the accretion process (§3.2). This means that the accretion is 
not effective for dwarf galaxies. Even if the efficiency of the accretion /5 acc , determined by 
the lifetimes of molecular clouds (D98), is as high as that in spiral galaxies, the accretion 
is not effective because of the low metallicity in the dwarf galaxies. Therefore, we cannot 
attribute the ineffectiveness of the dust accretion process to the lack of molecular clouds. 

Finally, we note that our model have satisfied one condition which any model must 
fulfill: The model has to explain the observation of nearby galaxies. Then, it becomes a 
matter of concern whether our model can explain the galaxies in the high-redshift Universe. 
For theoretical modeling of the cosmic dust mass, see, e.g., Edmunds & Phillipps (1997). 
Observationally, it is interesting that high-redshift galaxies found recently have evidences of 
dust extinction (Soifer et al. 1998; Armus et al. 1998). The number count of galaxies in the 
far- infrared and submillimeter wavelengths, where dust reprocesses stellar light, is another 
interesting theme concerning high-redshift dust (e.g., Takeuchi et al. 1999). 
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FIGURE CAPTION 

FIG. la — The relation between the dust-to-gas ratio (T> ) and the oxygen abundance (X Q ) 
for various /i n> o (the condensation efficiency of dust from oxygen atoms). The data points 
for dlrrs (filled squares) and BCDGs (open squares) are from LF98. The other parameters 
are fixed to /3 acc = 2/3 S n = 10, a = 1, and 5 — 1. The solid, dotted, and dashed lines 
represent different values of / in0 (0.1, 0.05, and 0.01, respectively). 

FIG. lb — The relation between /o = T>q/Xq (the fraction of oxygen in the dust phase) 
and Xo (oxygen abundance). The values of the parameters and the meanings of the data 
points are the same as Fig. la. 

FIG. 2a — The same as Fig. la but for the different parameter sets (/i n ,o — 0.05, a — 1, 
5=1 and various /9 acc and /?sn)- The solid, dotted, and dashed lines represent the cases of 
/3 acc = 2/5sn = 5, 10, 20, respectively. The long-dashed line shows the case of no accretion 
process onto the preexisting dust grains (/3 acc = and (3sn = 5). The data points are 
identical to Fig. la. 

FIG. 2b — The same as Fig. lb but for the parameter sets identical to Fig. 2a. The data 
points are identical to Fig. lb. 



